Cosmological evolution of general scalar fields in a brane-world cosmology 
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We study the cosmology of a general scalar field and barotropic fluid during the early stage of 
a brane-world where the Friedmann constraint is dominated by the square of the energy density. 
Assuming both the scalar field and fluid are confined to the brane, we find a range of behaviour 
depending on the form of the potential. Generalising an approach developed for a standard Fried- 
mann cosmology, in Q, we show that the potential dependence V((f)) can be described through a 
parameter A = — V' /(y/HV), where m 5 is the 5-dimensional Planck mass, H is the Hubble 

parameter and V' = ^x. For the case where A asymptotes to zero, we show that the solution 
exhibits stable inflationary behaviour. On the other hand if it approaches a finite constant, then 
V{4>) oc -^j. For A — *• oo asymptotically, we find examples where it does so both with and without 
oscillating. In the latter case, the barotropic fluid dominates the scalar filed asymptotically. Finally 
we point out an interesting duality which leads to identical evolution equations in the high energy 
p 2 dominated regime and the low energy p dominated regime. 

PACS numbers: pacs: 98.80.Cq 
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I. INTRODUCTION 

Scalar fields play very important roles in cosmology, es- 
pecially in the early stage of the universe. The Inflation- 
ary scenario relies on the potential energy of a scalar 
'inflaton' field to drive a period of early universe accelera- 
tion. It is thought that the early universe could well have 
also been characterised by a series of phase transitionSjin 
which topological defects could have been formed 0, U • 
Moreover in context of the string theory, the natural val- 
ues of the gauge and gravitational couplings in our 4d 
universe are explained by the dynamics of 'moduli' scalar 
fields. 0, |||. More recently, potentials of light scalar 
fields have been invoked as a possible source of the Dark 
Energy which appears to be dominating the dynamics of 
our Universe today. The success of this approach has of- 
ten been established by demonstrating the existence of 
attractor solutions by investigating the asymptotic dy- 
namics of the system being studied. For example in the 
context of usual 4D Einstein gravity, particular poten- 
tials have been discussed 0, IS III El EI 0; EL whereas 
others have adopted an approach which does not specify 
a particular potential QlEHEl- However, recent higher- 
dimensional unification scenarios based on brane-worlds 
suggest that the gravitational law could be different from 
Einstein's during the early stages of our Universe. In the 
brane-world scenario, our Universe is a four dimensional 
hypersurface embedded in higher dimensions. Standard- 
model particles are confined to the brane, while grav- 
ity can pro pag ate in the higher-dimensional bulk space 

EE E3, EE tlm HE H3, El 

Both the models based on Horava-Witten and Randall- 
Sundrum II (R-S II) are interesting because they lead 
to a new type of compactification with gravity. They 
also make concrete predictions of how conventional grav- 
ity will be modified at high-energy scales. Many au- 
thors have discussed the cosmology associated with these 



scenarios [23 H \X& M, S3, IH HJ. When we con- 
sider the homogeneous and isotropic cosmology based on 
these brane world models, the difference from conven- 
tional cosmology can be related to the appearance of two 
new terms in the Friedmann equation, i.e., the quadratic 
term of the energy-momentum and a dark radiation term 
[3(|0nilHl The allowed im pac t of dark radiation, is con- 
strained by nucleosynthesis |33| , with the most important 
change in the scalar field dynamics being due to the ap- 
pearance of the quadratic energy density term. 

In this paper, we study the dynamics of a scalar field 
which is confined to the brane in the context of a R-S 
II brane-world. As well as the scalar field, we include 
a barotropic fluid on the brane (which could be radia- 
tion or matter). Since the modification of the cosmic 
expansion law arising through the quadratic energy den- 
sity term becomes important during the early stage of 
the universe, we concentrate on this stage. In order to 
classify the asymptotic behaviour of the solutions, we fol- 
low the model independent approach as first proposed in 
P, 0| . Earlier work investigating the impact of these p 2 
corrections has concentrated on models with particular 
potentials [M IM IM SI B III III We intend to go 
beyond that in this paper, and in passing note that the 
analysis complements our previous discussion on scaling 
solutions in p 2 cosmologies pi) . 

The organization of the paper is as follows. In Sec. II, 
we review the key equations of motion and define vari- 
ables which will allow us to analyse the asymptotic be- 
haviour of this system and point out an intriguing dual- 
ity with the equations obtained in the context of scaling 
solutions in standard Friedman cosmologies 0, |nj . In 
Sec. Ill we obtain the class of attractor solutions which 
exist for constant A. We then extend the analysis to the 
cases where A evolves and find classes of solutions corre- 
sponding to inflationary regimes (Sec. IV); scalar field ki- 
netic energy dominated regimes (Sec. V), and oscillating 
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regimes (Sec. VI). We summarise our results in Section 
VII. 



II. EQUATIONS OF MOTION 

We analyze the dynamics of a scalar field in a Randall- 
Sundrum II (R-S II) brane- world [2^ , because the model 
is simple and concrete. However, we expect our key re- 
sults to also hold in other brane- world models, in which 
a quadratic term in the energy-momentum tensor gener- 
ically appears. In the R-S II model, even though the 
extra-dimension is not compactified, gravity is confined 
to the brane at low energy, resulting in Newtonian grav- 
ity in our world which is described by the intrinsic metric 
of the 4-dimensional brane spacetime. By use of Israel's 
junction condition and assuming Z 2 symmetry, the grav- 
itational equations on the 3-brane are given bv|30j 



(4) G _ 



n 



.(4) 
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in this paper, we must consider the term arising from 

For matter fields on the brane, we consider a scalar 
field <j> with a potential V(4>) as well as a barotropic fluid 
with an equation of state Pb = (7 — 1)ab, where 7 is an 
adiabatic index. (We will mainly be considering the case 
of radiation with 7 = 4/3, since we are interested in the 
early stage of the universe before nucleosynthesis.) We 
will also assume that the two fluids do not couple to each 
other explicitly. 

Although a 5-dimensional scalar field living in the 
bulk|43j may also appear in other brane-world scenar- 
ios, in this paper we only consider a 4-dimensional scalar 
field confined to the brane. Such a scalar field might 
originate in the condensation of fermionic matter fields 
confined to the brane. 

As the universe expands, the energy density decreases. 
This means that the quadratic term could have been very 
important in the early stages of the universe. Comparing 
the energy density terms in Eq. we find that the 
quadratic term dominates when 



(1) 



where ^Gn U is the Einstein tensor with respect to the 
intrinsic metric g^, is the 4-dimensional cosmo- 

logical constant, represents the energy- momentum 
tensor of matter fields confined to the brane and n^,, 
is the local correction term as a result of a brane em- 
bedded in the bulk. is a part of the 5-dimensional 
Weyl tensor and carries information about the bulk ge- 
ometry. n 2 = 8irG4 and k 2 = 87rG5 are 4-dimensional 
and 5-dimensional gravitational constants, respectively. 
In what follows, we use the 4-dimensional Planck mass 
7774 = k^ 1 = (2.4 x 10 18 GeV) and the 5-dimensional 

Planck mass 11J5 = k 5 2 ^ 3 . 771,5 is related to 7774 in terms 
of the brane tension a as am 2 = 67775, and could be much 
smaller than 7774. 

Assuming the Friedmann-Robertson- Walker spacetime 
in our brane world, we find the effective Friedmann equa- 
tions arising from Eq.Q arc 
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where a is a scale factor of the universe, H — a/a is its 
Hubble parameter, k is a curvature constant, P and p are 
the total pressure and energy density of matter fields, re- 
spectively. C is a constant related to "dark radiation" 
coming from |3ll l32|. In what follows, we consider 
only a flat Friedmann model (k — 0). We also assume 
that ( 4 ^A vanishes at least in the early stage of the uni- 
verse and C is negligible since it should be diluted by 
inflation and is strongly constrained at the time of Nu- 
cleosynthesis |23. (However, it must be noted that if we 
consider perturbations around the back ground solutions 



p> p c = 12777 5 6 /m 4 2 . 



(4) 



In this regime the expansion law of the Universe is 
clearly modified. For example, the expansion law in 
the radiation-dominant era becomes a cx i 1 / 4 compared 
to the conventional case, a cx t 1 / 2 . If such a period 
of p 2 domination did exist, then Nucleosynthesis would 
have provided us with a constraint on 7775, since it must 
have occured during the conventional-radiation domi- 
nated era in order to explain the abundances of the 
light elements. For example, assuming that the en- 
ergy density at a — a c is dominated by radiation as 
Pc ~ ( 7r2 /30)(7T ( 4 , where g is the number of degrees of 
freedom of light particles, the temperature of the uni- 
verse T c must be higher than that of nucleosynthesis, 
i.e. T c > Tjvs ~ lMeV. This constraint implies that 
7775 > 1.6 x 10 4 (3/100) 1 / 6 (T NS /1 MeV) 2 / 3 GeV. (For the 
brane tension a, this constraint is cr 1 / 4 > 2.0MeV). 

In the rest of the paper we will be concentrating on the 
evolution of the universe in the regime dominated by the 
p 2 contribution to the energy density. We will find there 
is a way of performing this analysis in an analogous man- 
ner to that already developed for the usual cosmologies 
dominated by the linear p term 0, 0, . In fact we will 
introduce a set of variables which will lead to exactly the 
same form of evolution equations for the two regimes, 
allowing us to relate the attractor solutions that exist 
in the two regimes through duality transformations. Of 
course, in a realistic cosmology, the p 2 dominated regime 
is eventually replaced by the linear regime. Any asymp- 
totic solution in the p 2 regime should be reached before 
the linear regime takes over at p ~ p c , and we confirm 
that is the case in this analysis. 

In the regime where p> p c , Eq. along with conser- 
vation of energy and momentum for each fluid component 
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on the brane gives 
1 



1 :r 



H = -j2^(pfl + 2^+V)(7p B + ^), (5) 



Pb = -Z~iHp B , 



dV 



(6) 
(7) 



which are subject to the Friedmann constraint given in 
Eq. © 



I' : - ^{PB + U 2 + V), 



(8) 



where dots denote derivatives with respect to time. The 
energy density and pressure of the homogeneous scalar 
field are given by p$ — 4> 2 /2 + V and — 4> 2 /2 - V, 
respectively. The effective adiabatic index of the scalar 
field at any time, 
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(9) 



is, in general, time or scale dependent. 

In order to analyze the stability of the solutions, we 
modify the variables, XclWj^CLW nrs t introduced in 
Ref. to investigate the evolution of a scalar field in a 
conventional Friedmann cosmology. Defining 



X 



Y = 



-3/2 1 
m 5 4> 

2\/3 VH' 

TTic V V 



V6 V7T 
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where a prime denotes a derivative with respect to the 
logarithm of a scale factor a, N = I n a, we see that these 
are related to those defined in Ref. by 



X 



Y k, H 



AxjLW ^CLW K4 V 2 



(12) 
(13) 



In terms of these new variables, the equations of mo- 
tion read, 

X' = -3X + X^Y 2 + ^X[2X 2 + 7(1- X 2 -Y 2 % 

(14) 

Y> = -\^-XY+ i -Y[2X 2 + 1 {l-X 2 -Y 2 % 

(15) 

A' = -V6X 2 {T-l)X + ^X[2X 2 +j{l-X 2 -Y 2 )}, 

(16) 

where generalising the expressions in Refs. Q and jg], we 
have defined 



A = -y/2 



m\ /2 dV/d(j) 
Vh V ' 



(17) 



r = v 



(PV/df_ 

{dV/d<t>) 2 



(18) 



Following Ref. 0, we see that the equation for H can 
also be written down as 



H' = -3H[2X 2 + 7 (1 - X 2 - Y 2 )}. 



(19) 



Since the energy density of the barotropic fluid satis- 
fies (pb > 0), the contribution of the scalar field to the 
total energy density, = p^j '(dm^H) = X 2 + Y 2 is 
bounded, < X 2 + Y 2 < 1. Hence the evolution of this 
system is completely described by trajectories within the 
unit disk. Moreover, since the system is symmetric un- 
der the reflection (X, Y) — > (X, —Y) and the evolution 
will not go beyond the Y = line (which corresponds to 
4> = 00 or H = 00), it is enough to discuss the upper half 
disk (Y > 0). 

There is an intriguing aspect to these evolution equa- 
tions. Given the definitions of X and Y above, Eqs. 1141 - 
I16fl are identical in form to those for Xqlw, ^CLW and 
Aclw = ~^~tt where A and Aclw are related through 



Aclw 



K4 2 
k 5 \ H 



(20) 



in P, 0| . This means that the attractor solutions 
for Xqlw etc, obtained in the standard Friedmann case 
where the Hubble parameter is determined by the linear 
energy density term, are the same attractor solutions in 
terms of our new variables X etc... but for the case of 
p 2 domination. Even though they do not correspond to 
the same potential or evolution of the Hubble parameter, 
they do correspond to the same parameters in terms of 
X, Y and A. This is intriguing, as it implies there is a 
duality between the high energy p 2 set of solutions and 
the low energy solutions corresponding to p domination. 



III. ATTRACTOR SOLUTIONS 

We want to now turn our attention to the evolution 
equations l|14|) - l|16fl . We are particularly interested in 
the case of scaling solutions, where X' = Y' = X' = 0. 
To start with we consider the case where A tends to a 
constant value asymptotically. Using Eq. (|15fl . for non- 
trivial X, Y and A, we see that the condition A' = 
in Eq. (|16fl corresponds to the constraint T = |. Inte- 
grating, Eq. (|18J) this in turn implies the form for the 
potential which leads to a constant A, namely 



V{<f) = 



(21) 



where p is a constant with the dimensions of mass. This 
result was first obtained using a different method in 
[4lj| . The equivalent potential which leads to scaling 
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solutions with constant Aclw i n the case of the stan- 
dard Friedmann cosmology is an exponential potential 
V(((>) — Vo exp(— Aclw k 4^) [ill - which is clear from 
Eqs. (|17|l and l|2(J[) . We can now summarise the attrac- 
tor solutions X c , Y c , for our scaling potential Eq. (|2T)l . 
recalling that they correspond to the same scaling solu- 
tions X C lw, Yclw as found in 0. Eqs. {HJ) and {I2J) 
admit five different critical solutions for X, Y with A 
constant. 

For these critical solutions, (X c = 1, Y c — 0),(X C = 
— 1, Y c = 0), (X c — 0, Y c — 0) are unstable (extreme) 
critical points. The other two solutions depend on the 
value of A. For A 2 > 37, we find the critical attractor 
values 



Yr 



2 A' 



3(2-7)7 



2A 2 



(22) 



1 x ? — 7 implying that both the 



which satisfies 7^ = x -± +Y 
barotropic fluid and scalar field energy densities scale 
with the same redshift. The fractional energy density 
stored in the scalar field is given by fl^ c = (37)/A 2 . 
If we consider the case of radiation for the barotropic 
fluid, then from Eqs. {TTJ, (0 and J22J) we obtain 
Q<j>c — (l/2)(^/ m 5) 3 |4l|- It is worth noting that because 
of the nature of the attractor solutions, for 37 < A 2 < 67, 
even if the scalar field dominates the universe, the cosmic 
expansion law becomes the same as the barotropic fluid 
dominated universe. 

If the asymptotic behavior of this system is character- 
ized by this scaling solution, then once the energy den- 
sity has decreased to p ~ p c , we recover Einstein gravity 
with the linear energy density term dominating over the 
quadratic term in Eqs. @ and It is well known that 
in this conventional cosmology, a scalar field with poten- 
tial V = /i 6 0~ 2 leads to a tracker solution and dominates 
the universe at late times 0, 0] ■ In order for this term 
to then provide the observed Dark Energy today, it re- 
quires /i to be fine tuned to a value O(GeV), meaning 
that f^ c during the period that the quadratic term dom- 
inates must be extremely small. 

On the other hand, if A 2 < 6, the corresponding at- 
tractor solution is 



X c 
Y. 



A 

7r 



1 - 



A 2 



(23) 



which satisfies 7^ = (1/3) A 2 and fl^ c — 1. In this case 
the scalar field energy density dominates the universe 
asymptotically. If the scalar field has 7^ < 7, then the 
solutions in Eqs. (|23|l are stable, and the redshift of the 
scalar field is slower than that of the barotropic fluid. 
It corresponds to power-law inflation and is stable as in 
the conventional cosmology driven by an exponential po- 
tential 0, 0| . Even though this solution is interesting 



because it is exact, it has the drawback of failing to stop 
inflating, so in a successful cosmology it would have to 
eventually end inflation and recover the conventional cos- 
mology. Furthermore, for such a large value of /i, this 
scalar field dominates the universe much earlier. How- 
ever, if 70 = (1/3) A 2 > 7, then the solutions in Eqs. Ij23(l 
are unstable and the scalar field ends up in the regime of 
the solution given in Eqs. (|22|) . 

Having considered the case of a constant A, we now 
turn our attention to the case where A(iV) is not con- 
stant. Then the critical values in Eqs. (|22H and (|23|l hold 
X 1 = Y' =0 only at single points, not over inter- 
vals of time. This means that the attractor solutions to 
Eqs. Ijl4(l and l|15|l are only valid as an asymptotic limit 
and in general X Cl Y c , A are time dependent. If X, Y 
do not oscillate, then since their values are constrained 
to \X\ < 1, |y| < 1, this implies that they will approach 
constant values asymptotically, given by the attractor so- 
lutions of Eqs. lfT4l) and (|T51) . and X\ Y' will vanish. 
Therefore, we can generalize the attractor solutions of 
X, Y given in Eqs. (1221) and (|23|l for more complicated 
potentials that have a nonconstant X(N). In the follow- 
ing sections, we will analyze other asymptotic solutions 
making use of this generalization. 



IV. POTENTIAL DOMINATED 
INFLATIONARY SOLUTIONS 



V2i 



3/2 



V 
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We now consider the case where A = — — ^ 
asymptotically. We might naively expect this happens 
as long as V — > faster than V as <j> — > 00, such as 
the case for inverse power law potentials of the form 
V(4>) — /i 4+ ™0~™, with n > 0. Indeed, in the conven- 
tional Friedmann cosmology in which Aclw = ~~kl^v 
plays the same role as our A, our expectation is realised 
in that Aclw — > 0, asymptotically [I]]. However in our 
case, if we substitute V((f>) — n A+n <f>~ n with n > into 
Eq. (|17|) . and use the fact that in the regime where the 
potential is dominating, H cx V, we obtain A cx 
From this, even with inverse power law potentials, only 
when 2 > n > 0, does A — > as 4> — * 00. 

In the following we will concentrate on models with 
potentials which lead to A — > asymptotically. If this is 
realized, we can eliminate the term proportional to A in 
Eq. I)14[l. and Eq. Q15[). After rewriting them in terms of 
H' with Eq. {EH, and since -6 < H'/H < for all values 
of X, Y, and 7, we obtain the following relations identical 
to those obtained in but recalling the interpretation 
for the variables X, Y and A differ in the two cases: 



= - 3 



X 

XL 
T 



1H'\ 
2F I<0 < 



1H' 



(24) 



This of course is to be expected given the duality we 
discussed earlier. As explained in [jj] these equations in- 
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dicate that X — > its minimum value, as Y —* 1, its max- 
imum value. This allows us to solve Eq. I|14|). Eq. (|15|l . 
and Eq. {El for -X", an d # m the region \X\ <C 1, 
|A| « 1 : 



X(N) 
Y(N) 



-3N 



ce -3 7 iV ' 
1 



ff(iV) = d(l 



-37W 



(25) 



where c and <i are integration constants. These approx- 
imate solutions are consistent with those in Eq. I|23[) in 
the limit A — > 0, whilst recalling that X and Y are vary- 
ing with time. They show that the scalar field potential 
dominates the energy density of the universe as the uni- 
verse inflates with almost constant Hubble parameter. 
In fig. 2] we demonstrate the validity of the solutions 
in Eq. I)25[l. as we show numerical solutions for a model 
with a barotropic fluid of radiation (7 = 4/3), and with 
potential V(<f>) = p 5 ^ 1 . 



(1a) 



(1b) 



(1c) 



(Id) 



^ Of 
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FIG. 1: 

the Friedmann equation. The universe contains radiation and 
a scalar field with V = p 5 ^ -1 . For simplicity, we choose 
fj, = 7715. The initial conditions are Xq — 0.2, Yq = 0.4. In (a) 
the Hubble parameter (solid curve) tends to a finite constant 
leading to inflation. In (b) we show Q^, quickly approaching 
unity. In (c) A evolves to zero asymptotically. In (d), the 
numerical solutions for Y (upper curve) and X (lower curve) 
are plotted. The solutions given by Eqs. 1251 are obtained 
asymptotically. 



This asymptotically inflating solution is stable, hence 
in order to recover standard cosmology at late times, a 
mechanism to end this period of inflation is necessary. 
Simply relying on the existence of the p term dominat- 
ing in the Friedmann equation will not be sufficient here, 



because as shown in 1], this potential also has an infla- 
tionary solution in that era. This can be readily seen 
from Eq. JSUJ, where Aclw = {^5 / K4) \/H/2 A — > as 
A — > for constant H . 



V. KINETIC-TERM DOMINATED SOLUTION 

Having considered the case A — * const (including 
zero), we now consider the case where asymptotically 

^2 m 3/2 , 

A = -7= — "yt - — > 00 smoothly without oscillations. 



Vh v 

The presence of the 1/y/H term in the definition of 
A, as opposed to the definition of Aclw V /V, and 
the fact that in most cosmologies H — ► asymptotically, 
means that a wider class of potentials satisfy A — * 00, 
than satisfy Aclw - * 00 • For example, as we show be- 
low, the inverse polynomial potentials, V{<j>) — p 4+n (f)~ n 
with n > 2 discussed earlier satisfy this condition. An- 
other concrete example is a model including an exponen- 
tial potential which satisfies A cx \j\fll — > 00. 

Although the specific potentials included in this limit 
may be different between the p 2 and p cosmologies, the 
invariance of the form of the equations of motion 114|) - 
(|15H when written in terms of X and Y implies that the 
solutions obtained in terms of X and Y in Ref. for the 
case of the conventional p dominated cosmology, apply 
to our case aswell. We therefore summarise the results 
presented there. 

Initially we expect X and Y to be order unity, which 
means for |A| ^> 1, the leading terms of Eq. 114|) and 
Eq. O are 



X' 



Y' = -\I~\XY. 



(26) 



Note for positive (negative) A, X — » 1(— 1) whereas for 
both signs, Y — > 0. Eventually as Y keeps decreasing, 
other terms in Eqs. (|14fl and l|I5|l become important. In 
particular if \X'\ > \X\Y 2 , then using Eq. I|19|) . the evo- 
lution of X is given by 



x' ( w 



(27) 



Now since —6 < H'/H < 0, X having reached a max- 
imum value turns over and like Y , heads off finally ap- 
proaching the values given by the solution of Eqs. I|22|l . 
X ->■ X c , Y -> Y c . For X c and Y c with |A| > 1, then 
1 3> \X C \ > Y c holds, and p$ will decrease faster than pb 
with — > 0. Even though X c , Y c are not critical (con- 
stant) points since A is not constant, we have verified nu- 
merically that the above asymptotic solutions are good 
approximations. In Fig. [2 typical examples of the nu- 
merical results for a model with potential V((j>) = p 7 4>~ 3 
are shown, in a universe containing radiation. Note in 
particular the behaviour of A — > 00 in (c) and how in (d) 
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X increases initially before decreasing towards zero, ap- 
proaching the solutions given in Eqs. As expected 
these solutions match closely those presented in Figure 
2 of Ref 0, for the conventional p dominated cosmol- 
ogy, but there are important differences. In particular 
the form of the potential involved (exponential in their 
case, inverse polynomial in ours), and most significantly 
the fact we are dealing with different cosmologies. 




FIG. 2: Evolution of the p 2 cosmology for a universe filled 
with radiation and a scalar field with V = p 7 cj>~ 3 - For simplic- 
ity, we choose p = ms. The initial conditions are Xo —0.2, 
Yo = 0.6. In (a) we show that the Hubble parameter normal- 
ized by its initial value (solid curve) tends to zero. In (b) we 
show quickly approaching 0. In (c) A evolves to infinity 
asymptotically. In (d), the numerical X (solid curve) and Y 
(dotted curve) solutions are plotted. As in Eqs. (I26L at the 
beginning, X increases and reaches its maximum value, be- 
fore turning over following Eq. 1271 as it (and Y) approaches 
asymptotically. 



Let us briefly digress to discuss the impact of this sort 
of solution. We have seen that as long as p < m 5 , then 
the asymptotic behaviour in the p 2 cosmology, is one 
where p 7 > p^. This will carry on until p ~ p c , after 
which we will enter a conventional radiation dominated 
universe. This is a nice feature as once in that regime we 
know from earlier studies that the inverse polynomial po- 
tentials being discussed here lead to tracking behaviour 
necessary to explain the dark energy today. The power 
of the p 2 regime here lies in the fact that it provides a dy- 
namical explanation of the large difference between the 
energy density of the dominant barotropic fluid and the 
scalar field 'initially' as we enter the conventional cos- 
mology |42|. 

If we decide to concentrate initially on the limit (f> — > 0, 
then this corresponds to A — > suggesting an unstable 



inflationary solution in this limit, for the inverse poly- 
nomial potentials V{4>) — p 4+n <fi~ n , with n > 2. This 
result provides one of a number of possibilities which 
have been investigated where the same scalar field po- 
tential is used to realize both a period of early inflation 
and late time dark energy domination, although there 
are strong constraints on their viability arising from the 
tendency to overproduce gravitational waves durin g re - 
heating. This is a form of quintessential inflation |46| . 
except that it is making use of the brane- world scenario, 
and has been invest igat ed recently by a number of au- 
thors m m m n m eh eh 

VI. OSCILLATING SOLUTIONS 

We now turn our attention to consider the case |A| — > 
oo with the scalar field <fi oscillating about the minimum 
of its potential asymptotically. In the context of realistic 
scenarios this is somewhat of a formal exercise in that 
we might well expect that the p 3> p c constraint would 
have been violated once this situation had been reached. 
However, we feel that it is worth investigating in its own 
right as it allows us to pursue the duality relation we 
have identified in the paper between the two regimes of 
p domination and p 2 domination. 

Given that <j) is oscillating about its minimum, we can 
without loss of generality, take it to be zero, and ex- 
pand V(4>) as a power series about the minimum. Keep- 
ing only the leading term, we have (1 / 'n) p A ~ n 4> n , where 
n = 2,4,6, •• •, even because of the boundedness of the 
potential. (For n — 4 we use V{4>) = (1/4)A0 4 where A 
is a dimensionless constant. ) 

Again because of the duality invariance, the equations 
match those in pj, so we follow their path to determine 
under which conditions, f2^: either dominates (goes to 
unity); oscillates around a finite constant value; or van- 
ishes asymptotically. In actual fact in |4l|, we have al- 
ready analyzed the cases for n = 2 and 4 obtaining 
analytic solutions by making assumptions about which 
components of the energy density dominate the universe, 
initially. In these cases, we have obtained approximate 
solutions by invoking the virial theorem, in which the 
relation between the time-averaged value of the kinetic 
term and the potential term of the energy density of the 
scalar field is given as 

(<j> 2 )=n(Vm. (28) 

From Eq. I|28|l . it is clear that the scalar field behaves 
as a pressureless perfect fluid (dust fluid) ((7^) = 1) for 
n — 2, while it behaves as a radiation fluid ((70) = 4/3) 
for n — 4. Therefore, if we consider the case of radiation 
for the barotropic fluid in the early universe, for the case 
n = 2, then the energy density in the scalar field will 
decrease slower than that of radiation, leading to f2^ — > 1. 
For the n = 4, case, both fluids evolve at the same rate, 
hence fi^ — » cie, a constant value determined by the self 
coupling parameter A in V(<fi). 
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In order to analyze the case for general n, and without 
specifying which component dominates the universe, we 
will introduce the total adiabatic index 7 tot . When it is 
the barotropic fluid that dominates, then 7 to t = 7, but 
when the scalar field dominates the universe, we take 
7tot = (7</>) since in this case 74, oscillates. Now, in terms 
of 7totj H(t) = 1/ (37toti), hence in the p 2 dominated era, 
Eq. (7J can be rewritten as, 



1 



4—nin—l 



+ -6 + M 4 -"0 

7tot* 



(29) 



(Recall for the case n = 4, the third term must be A</> 3 . In 

that case, if we introduce a new time variable t = VAt, 
the differential equation takes the same form as Eq. i|29|) 
with n = 4. Where appropriate we use this time variable 
for the case n = 4.) 

Even though we have already obtained the results for 
the case n = 2 in we first consider this case to con- 
firm our previous results. For V((f>) — (1 / '2) p 2 </> 2 , the 
solution of Eq. I|29|l is given in terms of J m and Y m , the 
Bessel functions of the first and second kind, respectively, 

4>{z) = pz m [a J m {z) + c 2 Y m (z)} , 

4>(z) = p 2 z m [c 1 J m ^ 1 {z) + c 2 Y m ^ 1 (z)], (30) 

where cu and c 2 are constants, m = (1/2) — (l/27 to t), 
z = pt 12]. From Eq. H30|) . X and Y are expressed as 



X = rz m+ ^[kJ m ^(z) + Y m ^(z)}, 
Y = rz m+ i[kJ m (z) + Y m (z)}, 



(31) 



where r = (l/2)(p/m 5 ) 3 / 2 v /j to tC 2 , k = ci/c 2 . A simple 
analytic expression can be obtained using the formula 
for the asymptotic limit of the Bessel functions J m ~ 
^(2/?rz)cos[z - tt(2to + l)/4], Y m ~ ^/(2/nz) sin[z - 
?r(2m + l)/4] for z > 1. The amplitude of X and Y in 
Eqs. (J3T} in the limit z > 1 goes as X ~ 7 ~ z m . A 
finite value of X, Y requires m = (1/2) — (l/27 to t) = 
(i.e. 7tot = 1)- Next, we analyze how the asymptotic 
value of 7tot are decided. In this case, regardless of the 
form of the barotropic fluid, in the limit 2> 1, (X 2 ) — 
(Y 2 ) = (r 2 /n)z 2m {k 2 + 1), which yields ( 70 ) = 1. This 
is completely consistent with the previous results relying 
on the virial theorem. We can confirm, therefore, if the 
barotropic fluid has 7 < 1, i.e. smaller than 7^, then 
the barotropic fluid decays slower than the scalar field, 
which yields 7tot < 1. In this case, from the asymptotic 
limit of the Bessel functions, fi^ = X 2 + Y 2 — > 0. On 
the other hand, if the barotropic fluid has 7 > 1, i.e. 
greater than 7^, then the barotropic fluid decays faster 
than the scalar field. In this case, 7 to t - ► (l4>) = 1, which 
leads to f2^ = X 2 + Y 2 — > 1 asymptotically. For the case 
7 = 1 = 7tot, ^4> asymptotes to some constant value 
between and 1. 

We have obtained an analytic solution of Eq. (|29|l for 
the case V(<fi) = (l/2)p 2 (f> 2 . This is clearly the simplest 



case since Eq. (|29|) is linear in 4> and its derivatives. For 
V{<p) = {l/n)^- n (p n with n > 2, Eq. ^ becomes non- 
linear in 4> and no simple analytic solution exists. How- 
ever, using the following ansatz, 



4> = [iz" 



C\ COS 



C2 sin 



(" + 2) 



, (32) 



we obtain a solution with the correct asymptotic behav- 
ior 0™/ 2 , cf> (x i" 1 / 2 , H cx t" 1 , and <fi oc neces- 
sary if we wish X, Y to have finite values other than 
and 1. The z(= /it) exponents in Eq. (|32() are such 
that Eq. (|29(l is satisfied up to leading orders i -3 / 2 and 
same condition also imposes the conditions 
7tot = 2n/(2 + n) and /3 2 = (27i/n + 2) 2 C2~ 2 [sin(/3z 2 /") + 
k cos(/3z 2 /")] n_2 . Notice that the value of 7tot is precisely 
the one we obtained earlier based on general arguments 
[cf. Eq. (|28p]. Of course, the ansatz in Eq. (|32|l is not a 
complete solution to Eq. I|29|) as can be seen in the fact 
that (3 is not a true constant (except for the case n = 2). 
For n > 2 we must use the average value of (3 evaluated 
in the asymptotic region. In terms of Eq. H32[). X and Y 
take the following expressions: 



X = X 
Y = Y a 



30s ( pz 2 ™ J — k sin I f3z 2 ™ J , 
m(/3z^^) +kcos((3z^^\ ,(33) 



W ith Xg = ^^(^/ m 5) 3/2c 2 aild Y = 

V / W^(/'/™5) 3/2 cf. 

From Eq. |j33I) we obtain (Y 2 )/(X 2 ) = 2/n, asymptot- 
ically, as in Eq. (|28|l . which is consistent with the virial 
theorem. We have, therefore, (7^) = 2n/(2 + n) and 
7tot = (70), i- e -, the ansatz in Eq. I|32l) is a solution to 
Eq. (|29J) only when the dominant energy density redshifts 
as fast as the scalar field. This is, of course, no surprise 
since we imposed on the ansatz Eq. (|32|) . the limit X, Y 
->■ cte (7^ or 1). 

Even though the above discussion strictly relies on 
the time-averaged value, we have confirmed that it pro- 
vides a good approximation to the true numerical re- 
sults. In fig. [31 typical examples of the numerical results 
for models with V((j>) = (l/2)p 2 4> 2 (in (a) and (b)) and 
V{<j>) = (1/4) \<j) 4 (in (c) and (d)) are plotted. As for the 
barotropic fluid, we use radiation (7 = 4/3). For both 
potentials, we draw the time evolution of and Ht. 

To conclude this part of the analysis, if the ini- 
tially dominant energy density component has a larger 
(smaller) adiabatic index than (7^) = 2n/(2 + n), then 
Q<f, will approach 1 (0). For example, if we consider radi- 
ation as the barotropic fluid, then for n > 4, the energy 
density of the scalar field will decrease faster than ra- 
diation, which guarantees the recovery of the standard 
cosmology at an epoch when the energy density is lower 
than p c . 

Even though the form of these relations are the same 
as in a conventional cosmology, the oscillating behavior 
found in the regime dominated by the quadratic term 
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(3b) 



A 



-Jh v ~ p 



70 



V{<f>) 



/ 



c = cte (> 
c = cte (< 
oo (nooscil.) 
oo (oscil.) 



37 
2c 5 " 

1 




cte 
1 
1 



7 

3 

7 



2+n 
2n 

2+n 
2n 



(>7) 
(=7) 
«7) 



/' 



6 0" 2 

n>2 
n > 0, even 

\ 2 > n > 



1 1 1 



(3c) 



15 2 2 5 24 




TABLE I: The asymptotic behavior of the contribution of the 
scalar field to the total energy density (£1$) and the effective 
equation of state for the scalar field (7^) for different limiting 
cases of X(4>) in a p 2 dominated cosmology. In the last column 
we give an example of the potential V((f>) which satisfies this 
limit. Note how the duality transformation implies solutions 
match those in pj for the variables A, X and Y, although 
the definitions of those variables differ. In particular, the 
corresponding potentials in the last column differ for the two 



FIG. 3: Evolution of the universe filled with radiation and 
a scalar field with potential V = (1/2) p 2 (j> 2 ((a), (b)) and 
V — (l/4)\(j) 4 ((c), (d)) in the p 2 dominating era. For sim- 
plicity, we choose p = 7715 and A = 1. The initial conditions 
are Xo = 0.1, Yq = 0.5. We confirm, as expected, that the 
asymptotic behavior of the oscillating scalar field is character- 
ized by (70) = 2n/(2 + n). In (a) we show the time evolution 
of flp with a (f> 2 potential (n = 2). Since the oscillating scalar 
field behaves like dust for n = 2, it gradually dominates the 
radiation and fi^ oscillates, asymptoting to unity. Similar 
behavior is seen in Ht in (b), where it approaches to 1/3 
corresponding to the cosmic expansion law of a dust domi- 
nated universe in the p 2 dominating era. In (c) we show the 
time evolution of O.^, for the (j) 4 potential, (n = 4). Since the 
oscillating field behaves like radiation here, after a few oscil- 
lations, it approaches a constant value. In (d), for n = 4, 
Ht — > i , corresponding to the cosmic expansion law of a dust 
dominated universe in the p 2 dominated era. 



has a distinct advantage in the preheating phase [44( if 
this scalar field couples to other matter fields. Since the 
cosmic expansion at late times becomes slower and the 
amplitude of (f> decreases more slowly than in the conven- 
tional case, there may well be sufficient particle produc- 
tion generated as the field decays. 



VII. SUMMARY 

In this paper, we have studied the novel dynamics of 
a scalar field plus barotropic fluid, in a brane-world sce- 
nario dominated by the p 2 term in the effective four- 
dimensional Friedmann constraint. As a concrete exam- 
ple, we have adopted the Randall- Sundrum II model and 



assumed that the scalar field is confined to our four di- 
mensional spacetime. 

Our approach has allowed us to deal with general 
classes of potentials, and complements an earlier inves- 
tigation of a similar system but for particular potentials 
in |4l|. Perhaps the most important result we have ob- 
tained, can be seen in the defining equations HlUfl - l|19|) m 
which we introduce a new set of variables to analyse the 
evolution equations in a model independent manner. The 
crucial point is that the equations for X and Y, given by 
(|14fl and H15fl are identical to those derived in the conven- 
tional four dimensional cosmology where the Friedmann 
equation is driven by the energy density p of the scalar 
field and barotropic fluid, Q, [ill IT^ as opposed to the 
p 2 driving term in our case. Although the precise def- 
initions of X and Y differ, the fact that they obey the 
same evolution equation allows us to immediately write 
down and understand the form of the scaling solutions. 
There is a simple map which allows us to relate the cos- 
mological solutions in each regime - a duality between 
the parameters X, Y in the two cases. 

We have summarized our results in Table I, in an anal- 
ogous manner to that presented in 0] for the equiva- 
lent conventional p dominated cosmology. In particular 
we showed that all the model dependence is given by 

\/2m 3/2 v' 

A = — and the adiabatic index of the barotropic 

fluid 7. Scalar potentials that do not require introducing 
nonzero minima are classified into one of three different 
limiting cases by the asymptotic behavior of A: A goes to 
a finite constant, zero, or infinity. In the first case, 
approaches a finite constant (different from one or zero) 
depending on the value of A. It is worth noting that in the 
conventional cosmology this happens in the model with 
an exponential potential, V(<fi) — /x 4 exp (— Aclw k 4^0j 
while in the brane-world cosmology, it happens with an 
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inverse square potential model, V{<j>) = p 6 <fi~ 2 - a reflec- 
tion of the duality relating X, Y and ^clWi ^clw- 

In the second case A — > 0, we obtained X — > 0, Y — > I 
corresponding to inflation with an almost constant Hub- 
ble parameter. As a concrete example, we considered 
the model with V(<f>) — /i 4+ ™</>~™ . For this potential, in 
the conventional cosmology, the inflationary solution is 
an attractor for any n > 0, but in the presence of the 
quadratic density term, it occurs only for 2 > n > 0. 

In the final case A — > oo. Following the procedure in [l| 
we investigated two different possibilities depending on 
whether A oscillated or not. If it didn't, X, Y, fl^ — ► 
asymptotically. We showed the model with an exponen- 
tial potential belongs to this class in the presence of the 
quadratic energy density term, as do models with inverse 
power law potentials with n > 2. As pointed out in |42j . 
this provides a new feature for the quintessence scenario 
in that it allows us to explain how the energy densities 
in the radiation and scalar field could be different as we 
enter the usual p dominated Friedmann cosmology era. 
On the other hand, if A does oscillate, the adiabatic index 
of the scalar field is given as (7^) = 2n/(2 + n), where 
n is the power of the leading term in the scalar poten- 
tial. The asymptotic behavior of the universe is decided 
by whether (7^) is larger than 7 or not. For (7^) > 7, 
Q,<p — ► 0, while for (7^) < 7, Q</, — > 1 and for (7^)7, 
Cl<p — > cte. 

Although, we have concentrated on investigating the 
cosmology in the p 2 dominated regime, any realistic cos- 
mology also has to allow for the fact that the universe 
is dominated today by the standard p term in the Fried- 



mann equation. The next thing to do is to match the 
two regimes together. This has been investigated by 
a number of authors [MIMIMIIlHIIallSIIll 
mainly using numerical simulations involving particular 
potentials. Such quintessential inflation models are being 
tightly constrained by recent CMBR data as they tend 
to generate too lar ge a tensor contribution to the CMBR 
power spectra [35|, |40( ■ What we have developed in this 
paper is an alternative, possibly powerful approach which 
allows us to make use of the duality that exists between 
the two regimes (p 2 and p dominated). We are currently 
investigating using this duality to determine consistent 
cosmologies involving an evolution from one regime into 
the other in a model independent manner, but based on 
the idea of using the definition of A and 7 as the key 
ingredients in the analysis. 
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